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ABSTRACT 

Radio loud active galactic nuclei (AGN) are on average 1000 times brighter in the radio band compared to 
radio quiet AGN. We investigate whether this radio loud/quiet dichotomy can be due to differences in the spin 
of the central black holes that power the radio-emitting jets. Using general relativistic magnetohydrodynamic 
simulations, we construct steady state axisymmetric numerical models for a wide range of black hole spins 
(dimensionless spin parameter 0.1 < a < 0.9999) and a variety of jet geometries. We assume that the total 
magnetic flux through the black hole horizon at radius rn(fl) is held constant. If the black hole is surrounded 
by a thin accretion disk, we find that the total black hole power output depends approximately quadratically 
on the angular frequency of the hole, P oc 0^ oc (a/rn) 2 . We conclude that, in this scenario, differences in the 
black hole spin can produce power variations of only a few tens at most. However, if the disk is thick such that 
the jet subtends a narrow solid angle around the polar axis, then the power dependence becomes much steeper, 
P oc or even oc fl^. Power variations of 1000 are then possible for realistic black hole spin distributions. We 
derive an analytic solution that accurately reproduces the steeper scaling of jet power with ilu, and we provide 
a numerical fitting formula that reproduces all our simulation results. We discuss other physical effects that 
might contribute to the observed radio loud/quiet dichotomy of AGN. 

Subject headings: relativity — MHD — galaxies: jets — accretion, accretion disks — black hole physics — 
galaxies: nuclei — quasars: general 



1. INTRODUCTION 

The first active galactic nuclei (AGN) were discovered 
\ through radio emission associated with their relativistic jets. 
However, it soon became clear that not all AGN 3 have pow- 
erful radio jets; in fact, only about 10% of quasars do. The 
evidence for a dichotomy between radio loud and radio quiet 
AGN has become stronger over the years, culminating in 
the im pressive demonstration by ISikora. Stawarz. & Lasotal 
I (120071) that two very distinct populations of AGN are clearly 
visible when radio luminosities Lr of AGN are plotted against 
optical luminosities Lb- For a given value of Lb, these authors 
show that Lr of radio loud AGN is ~ 10 3 — 10 4 times greater 
than that of radio quiet AGN. Also, the two populations fol- 
low two well-separated tracks on the plot. 

The origin of the radio loud/quiet dichotomy has been 
much discussed in the literature. At Eddington ratios A = 
iboi/^Edd ~ 10Lfl/LEdd > 0.01, where Lboi is the bolometric 
luminosity of the AGN and Lem is its Eddi ngton luminos- 
ity, a likely expla nation for the dichotomy dHo et al.l 120001: 
ISikora et alj|2007l) is that these systems accrete via a standard 
thin accretion disk. Jet production is then expected to be in- 
termitten t, as found to be th e case in black hole (BH) X-ray 
binaries dFender et al.ll2004h . However, the existence of two 
distinct populations for A < 0.01 is harder to explain. Even 
at these low luminosities, the radio loudness parameter R = 
L5 ghz/ Lb of the radio loud population is at least a factor of 
10 3 times larger than that of the radio quiet population. How- 
ever, at low values of A, BH X-ray binaries typically are in 
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a hard spectral state as sociated with an advection-do minated 
accretion flow (ADAF, iNaravan & McClintockll2008l). and in 
this s tate, all BH X-ray binaries are radio loud dFender et al.l 
120041) . Why then do AGN with similar values of A have a 
radio loud/quiet dichotomy? 

One possible explanation is that radio loud objects are 
driven by a central BH with a large spin which pro- 
duces a jet by the Bla ndford-Zn ajek ( BZ) mechanism 
dBlandford & ZnaieklfT977l hereafter IBZ77I). This is referred 
to as t he spin paradigm dBlandford! 1 1990b I Wilson & Colberfl 
11995b iBlandfordl 119991) . which is in contrast to the accre- 
tion paradigm which states that the BH mass and mass ac- 
cretio n rate determine the jet power dBlandford & Reesll 19741 
119921) . These d ifferent paradigms plausibly operat e together 
to some degree dBegelman et al.lll984tlMeierJl2002l) . The spin 
paradigm has been invoked to explain the observed correla- 
tion between jet and accretion power in elliptical galaxies 
d Allen et al.1 120061) b y combining an ADAF accretion model 
with the BZ effect dNemmen et al.l 12007b iBenson & Babull 
120091) . In terms of the dimensionless spin parameter a = 
J/GM 2 , where M and J are the mass and angular momen- 
tum of the BH, it is found that one requires a > 0.9 to explain 
the correlation. The dichotomy in the power and spatial dis- 
tribution of emission in Fanaroff-Riley classes 1 and 2 (FR 1 
and FR 2) radio galaxies m ay also be explained by the spin 
paradigm dBaum et al.ll 19951) . 

The radio loud/quiet dichotomy in AGN dKellermann et al.l 
119891: iModerski et al.l 1 1998b llvezic et al.l 120041) has been ex- 
plained in terms of an in situ trigger for relativistic jets 
dMeier et al.l 1 19971) . It could also be explained by the dif- 
ferences in_the_eyolutionary stages at which we observe the 
objects dBlundellll2008l) coupled with the episodic activity of 
the AGN evidenced by the change or precession of jet orienta- 
tion dSteenbrugge & Blundellll2008b ISteenbrugge et al.ll2008l 
120101) . However, another possibility is that the two AGN pop- 
ulations have different merger and accretion histories which 
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lead to different BH spins. Recent observations show that, 
for A < 0.01, all radio loud AGN reside in ellipti cal galaxies, 
where a s radio quiet AGN live most l y in sp irals dSikora et al.l 
120071) . IVolonteri. Sikora. & Lasotal (120071) explored a num- 
ber of scenarios for the formation of ellipticals and spirals 
and showed that it is plausible for the nuclear BHs in spirals 
to have lower spins than those in ellipticals. This suggests 
that the spin paradigm may explain the radio loud/quiet di- 
chotomy. 

The main problem with this explanation is that the differ- 
ence in radio loudness between the radio loud and r adio quiet 
populations is a factor of 10 3 (ISikora et al.l 120071) . This is 
a strikingly large difference. According to the accretion or 
spin paradigms, relativistic jets are produced by magnetic out- 
flows from either the inner region of the disk or the spinning 
BH. The power in the disk outflow is expected to be propor- 
tional to the disk luminosity, which leaves no room for a ra- 
dio loud/quiet dichotomy, so we will ignore this possibility 4 . 
The power from the BH does depend on the spin parameter, 
and we will focus on this 5 . However, the analytical model 
of IBZ77I indicates that the jet power P varies only as a 2 for 
fixed magnetic flux threading the horizon. If such a weak 
variation is to produce a difference of 10 3 in radio power 6 , 
we need a in the two populations to differ by a factor ~ 30, 
which does not seem plausible given likely mer ger histories 
dHughes & Blandfordl20"0llGarrimie et al.l2004l) . More plau- 
sible is a factor ~ 3 differ ence in the median val ues of a in the 
two populations (e.g., see IVolonteri et al.ll2007l) . but this will 
produce only a factor ^10 difference in jet power, not 10 3 . 

The IBZ77I scaling for power, P oc a 2 , was derived in the 
limit a<l, for a razor-thin disk, assuming the magnetic flux 
threading the BH is independent of a. Recen t analytical and 
numerical work by Tanabe & Nagataki (2008) show that, for a 
BH threaded by a split monopole magnetic field, P increases 
as a 4 at large values of a when higher-order corrections are 
included. However, the analytical model worked out by these 
aut hors on ly gives a factor of two increase in power above 
the IBZ77I result even at a = 1. Their numerical simulations 
achieve a slig htly ste eper scaling, but still only a factor of 
four above the IBZ77I result at a = 1 . In any case, their work 
hints that a much steeper dependence of power on a occurs 
as a — > 1 . Are there other effects that can introduce an even 
steeper dependence on al 

General relativ istic magne t ohydro dynamic simulations of 
accretion disks bv lMcKinnevi d2005bl) showed that for a > 0.5 
the jet power varies as steeply as the fourth power of the BH 
angular rotation rate, i.e., P oc £1^, where i7e oc a/re and re 
is the radius of the horizon. Compared to the scaling P oc a 4 , 
the scaling P oc fl^ introduces an additional factor of 16 due 
to the division by rH, since r» decreases from 2M to M as 
a varies from to 1. iMcKinn ey (2005b) also finds that the 
power output of the entire BH has a shallower dependence on 

4 The BH spin also drives power into the disk causing a more powerful 
disk outflow, but this still requires BH spin to introduce a dichotomy. 

5 There are arguments to suggest that the luminosity of the disk 
outflow should be greater than that from the central spinning BH 
I Ghosh & Abramowicz 1997; Livio et al. 1999). However these arguments 
assume low values of turbulent viscosity and weak magnetic fields near the 
BH, and also do not account for the effects of th e general relativist ic plunging 
region (see, e.g., McKinnev & Naravan 20073 for a discussion). IMcKinnevI 
42005b!) finds that the lu minosities of the jet and the disk are similar (how- 
ever, see Nagataki 2009). For the purposes of this paper, we ignore the disk 
wind. 

6 We assume that the radio luminosity is proportional to the jet power. 



flu compared to the power output of the jet, which subtends 
a small solid angle above the disk and corona 7 . This sug- 
gests that changes in the solid angle subtended by the jet (via 
changes in the disk thickness) could change the steepness of 
the power output as a function of a. 

Since the scaling of P with BH spin is important for jet 
studies, we have explored this issue in detail using general 
relativistic magnetohydrodynamic numerical simulations. We 
consider a variety of geometries for the shape of the jet to see 
if we can come up with any scenario in which jet power could 
change by a large factor for a modest variation in a. We show 
that the most favorable scenario is a BH surrounded by a thick 
accretion flow with an angular thickness H/R ~ 1. We show 
that in this case the power output into a polar jet has a steep 
dependence on the spin, P oc il^, and that the scaling steepens 
to P oc il^ for e ven thicker disks. H ence, we confirm the basic 
result found by McKinnevl (l2005bl) of a steep dependence of 
jet power on a at high latitudes above the disk. We suggest 
that this strong dependence may explain the radio loud/quiet 
dichotomy in AGN. 

Our numerical setup is described in §[ 2] T he results for BHs 
with razor-thin disks are prese nted in 33.11 and for jets from 
BHs with thick disks in 33.21 We discuss the results in the 
context of the AGN radio loud/quiet dichotomy in §@] and 
conclude in §0 We work with Heaviside-Lorentzian units and 
set c = G = 1 . 

2. NUMERICAL SETUP 

It is known that a highly magnetized relativistic jet does 
not easily self-collimate. For instance, the equilibrium field 
configuration around an isolated spinning BH threaded by 
a magnetic field (sourced by ext ernal currents) takes the 
form of a split monopole (IBZ77I) . Only extremely close 
to the polar axis is any evid ence of self-collimation present 
dTchekhovskov et ai1l2009al) . Therefore, in order to produce 
a jet which collimates most of the energy output from the BH, 
it is necessary to introduce an external confining medium. The 
confining agent may be the gas in an accretion disk, a corona, 
or a wind emerging from the inner regions of the disk. Ideally, 
one should numerically simulate both the jet and the confin- 
ing medium, but this is numerically very ch allenging. Instead, 
we follow the more usual approa ch (e.g., iKomissarov et al.l 
120071: iTchekhovskoy et afll2009bl) of introducing a rigid ax- 
isymmetric wall with a prescribed shape and requiring the jet 
to lie inside the wall. The shape of the wall is set by two 
parameters: 

1 . An index v which sets the asymptotic poloidal field line 
shape, as described below. This parameter ranges from 
v = 0, which corresponds to a monopole field geometry, 
to v = 1, which corresponds to a paraboloidal jet. In a 
real system, v would be set by the radial pressure profile 
of the confining medium. Plausible values are in the 
range v ~ 0.5-1 dTchekhovskov et alj|2008l) . 

2. A transition radius ro which is defined such that for 
r < ro the jet is monopolar and for r > ro it switches 

7 McKinney 's models have an accretion disk with a disk+corona+wind of 
angular extent H/R ~ 0.6. For a > 0.5 the power per unit mass accretion rate 
scales as <x for the polar jet and oc for the entire horizon. However, in 
these models, the mass accretion rate through the BH horizon per unit fixed 
mass accretion rate at large radius scales a s 1/Qh for a > 0.5, beca use of the 
ejection of a massive wind (as also seen in Hawley & Krolik 12006)). Hence, 
expressed in terms of M at large radius, the polar jet power scales oc £1^ and 
the power output from the entire horizon scales oc . 




Figure 1. Meridional cuts through a selection of numerical jets in the force-free approximation, i.e., neglecting plasma inertia. The color-coded radial four- 
velocity u r (as measured in Boyer-Lindquist coordinates, see legend) is shown overlaid with poloidal field lines (thin black lines which correspond to if 1 / 2 = 
0,0.2, . . . , 1). The jets are confined by collimating walls (thick black lines). The left-most column of panels shows models with a low BH spin (a = 0.1), the 
right-most column models with a fast spin (a = 0.9999), and the middle column models with an in-between spin (a = 0.9). The top row of panels shows the most 
collimated paraboloidal models (v = 1 , rg = M), the bottom row non-collimating monopolar models (v = 0), and the middle row models that are monopole-like 
until r ~ ro = 10M, beyond which they become paraboloidal (v = 1 , rg = 10M). 



to the shape prescribed by v. The parameter ro allows 
us to consider situations in which confinement operates 
only beyond a certain distance from the BH. 

In terms of these two parameters, the wall has the following 
shape in polar (r, 9) coordinates in the Boyer-Lindquist frame: 



1 - cos 9 = 



r+r Q 



r H =M 



l+(l-a z ) 



2x1/2 



(1) 



where re is the radius of the BH horizon. For r<ro, (1 — 
cosf?) <C 1, so 9 « 7r/2, i.e., the wall lies along the equatorial 



plane, as for a split monopole. For r (rn,ro), 9 oc r 



-I//2 

which corresponds to a generalized paraboloid. Note that, in 
all these models, the wall meets the horizon at the equator. In 
effect, this means we assume a razor-thin disk which subtends 
zero solid angle at the BH. In §3.2 we discuss the case of 
geometrically thick disks. 

Having picked the shape of the wall, we choose the poloidal 
flux function of the initial magnetic field configuration to be 



<Jr = 



r+r 



(1-COS0). 



(2) 



Note that is conserved along each field line, and $ = 
27r<I'(r, 9) is the poloidal magnetic flux through a toroidal ring 
at (r, 9). By construction, equation (O corresponds to a total 
flux of <l> tot = 2ir in the jet. In this model $ tot does not depend 
upon spin, so the amount of magnetic flux threading the BH is 



fixed for different spins. The outermost field line, defined by 
\E f = 1, follows the shape of the wall. This particular field line 
is always located at the wall because of our boundary con- 
ditions. Interior field lines, however, are free to move once 
the simulation begins and do experience minor shifts in the 
poloidal direction. The initial magnetic field has no toroidal 
component. 

There are no known exact solutions for the magnetosphere 
of a spinning BH. However, the poloidal field configurations 
given in equation (fJI are sufficiently close to the true solutions 
that their initial relaxation when the simulation starts is rather 
mild. Only two linearly independent analytic solutions have 
been obtained for a non-spinning BH: one corresponds to a 
monopolar field geometry and is given by equ ation (f2| with 
v = 0, while the other is the following IBZ77I solution for a 
paraboloidal field, 



^_ (r/rH-l)(l-cos0)-(l+cos0)log(l+cos0) 
2 log 2 



(3) 



Note that the split-monopole solution applies to the entire 
space exterior to the horizon, whereas the paraboloidal so- 
lution only applies to the fiel d lines attached to the BH (see 
iMcKinnev & NaravarJ2007bl for a numerical paraboloidal so- 
lution that applies to the whole space). Any linear combina- 
tion of the mon opolar and par aboloi dal sol utions is also a so- 
lution (see, e.g.. lBesk in 2009). The IBZ77I paraboloidal solu- 
tion (O is very similar to the approximate solution (O for the 




Figure 2. Meridional cuts through a selection of numerical magnetohydrodynamic, i.e., mass-loaded, jets. These may be directly compared to the force-free 
(i.e., neglecting plasma inertia) jet models shown in Figure[T] the jet geometries are the same, but the jets here include mass-loading. While the field line shape 
changes little from Figure [JJ the introduction of mass-loading does lead to clear differences in the velocity, (i) The inflow near the BH is faster in mass-loaded 
jets (u r is more negative in blue regions) due to the gravitational pull of the BH. (ii) The outflow from the BH in monopolar jets (the bottom set of panels) is 
slower than in the corresponding force-free simulations. Nevertheless, the jet powers are nearly the same in the two sets of models. 



case v = 1, and this is our reason for focusing on the simpler 
model (|2), with v varying over the range to 1. For com- 
pleteness, we have also run simulations using the |BZ77| field 
geometry (0 to initialize the calculations (along with the ap- 
propriate choice of the wall shape, obtained by setting \E' = 1 
in this equation). 

We performed the simulations us ing the general rel- 
ativistic (GR) MHD code HARM (iGammie et all 120031: 
McKinnev & Gammid200llMcKinnevl2006bHatlNoble et al.l 
2001 " using Kerr-Schild coordinates in the Kerr met- 
ric; the code includes a num b er of recent improv e ments 
(iMignone & McKinnevI 120071: iTchekhovskoy et al.l I2007L 
120081 I2009al)~ Most of the simulations were done in the 
force-free approximation, which assumes that the plasma is 
infinitely magnetized and has negligible inertia. Within this 
approximation, the problem is fully defined by specifying just 
the BH spin and the shape of the wall. Figure Q] shows results 
from several representative simulations. 

Real relativistic jets are of course not perfectly force-free; 
in fact, they are expected to deviate substantially from this ap- 
proximation at large distances from the BH. However, all rela- 
tivistic MHD jets that have 7 ^> 1 are highly magnetized near 
the BH, and here the y are e xpected to be well represented by 
force-free solutions dBZ77h . Moreover, the power that a rel- 
ativistic jet carries is determined entirely by the initial force- 
free zone. Therefore, we expect numerical results on the jet 
power from a force-free simulation to agree very well with the 
power for an MHD jet with inertia. 



To verify this expectation, we have repeated some of our 
force-free simulations in the MHD limit, in which the jet is 
mass-loaded with a finite amount of plasma (details given be- 
low). Figure [2] shows some results. As expected, we find 
that the asymptotic Lorentz factor 7 of the jet does depend 
on mass-loading: a force- free jet accelerates without limit 
dTchekhovskov et ai1l2008l) . whereas an MHD jet asymptotes 
to a finite 7 which is determined by the initial magnetization 
of the jet. However, we find that the jet power, the primary 
quantity of interest to us in this paper, is not sensitive to mass- 
loading so long as the jet is relativistic, i.e., so long as the jet 
is initially force-free near the BH. 

MHD jets are more complicated and require more param- 
eters to be specified compared to force-free jets. In particu- 
lar, the results depend on the details of mass-loading at the 
base of the jet. Highly magnetized jets accelerate because the 
magnetic energy flux, which dominates the energy budget at 
the base of the jet, is converted to kinetic energy flux of the 
plasma as the jet flows out. The ratio of magnetic to kinetic 
energy flux at the base of the jet determines the asymptotic 
Lorentz fac tor (Bege lman & Lil l994tlKomissarov et alj2009t 
ITchekhovskoy et al.ll2009allbh . Observations suggest a char- 
acteristic value for t he Lorentz factor of AGN jets 7agn ~ 25 
dJorstad et al.ll2005h . We choose the following simple pre- 
scription for the mass-loading of our numerical MHD jets 
to roughly match this value. We impose a floor on the co- 
moving rest-mass density of the jet pfl oor such that when- 
ever the density falls below this value, we add mass in the 
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Figure 3. Jet power output of various models as a function of dimensionless BH horizon spin frequency ojh = S1h/^h(o =1) (equation[5j. Different models 
are shown with different symbols (see legend). For reference, the following values of BH spin, a = {0.2,0.5,0.9, 1}, correspond to approximately the following 
values of u)\\: {0.1,0.27,0.63,1}. [Upper panel (a)]: Logarithm of jet power. [Lower panel (b)]: Fractional deviation of the numerical jet power from the BZ6 
formula (5). [Insets (ia) and (ib)] Blow-ups of panels (a) and (b) for hig h values of BH spin. The values of the BH spin a are shown as labels next to datapoints. 
The three colored stripes (see legend) correspond to the original BZ77 scaling (equation|4), the BZ2 scaling (equation|7J, and the BZ6 scaling (equation|9). The 
second-order BZ2 scaling Q follows the numerical results very well except for values of a close to unity. The BZ6 scaling {9) matches the numerical data well 
at all values of a. For each a, the monopolar model (y = 0) has the maximum power and the paraboloidal model has the least power. However, the difference is 
only ~ 20% 



co-moving frame of the plasma. This simple floor model is 
a convenient way of numerically representing more compli- 
cated (and poorly understood) processes that are responsible 
for the mass-loading of jets in AGN. The value of pfl oor is se- 
lected such that the rest mass energy density pft oor c 2 is equal to 
a fraction 1 / 7agn of the co-moving magnetic energy density 
Cm- PfloorC = Cm /7AGN- Thus, our floor model ensures that the 
ratio of e m to Pfl oor c 2 does not exceed 7agn> thereby making 
sure that the maximum Lorentz factor of our jets is close to the 
required value. We note that while our procedure is Lorentz 
invariant, it might not correspond to a physical process that 
oper ates in AGN, e. g., photon annihilation from the accretion 
disk (lPhinnevll983l) . 



The code uses internal coordinates (x\,X2) that are uni- 
formly sampled with 512 x 128 grid cells. The internal co- 
ordinates are mapped to the physical coordinates (r, 9) via 

1 /2 

r/M = /?o + exp(xi) and X2 = sign(\I/) ' . The computa- 
tional domain extends radially from the inner boundary at 
r m = 0.6M+0.4r H to the outer boundary at r out = 10 4 M. We 
apply absorbing (outflow) boundary conditions at each of 
these boundaries. In the (^-direction, the computational do- 
main extends from the polar axis at = 0, where we use the 
usual polar boundary conditions, to the jet boundary X2 = 1, 
where we place the wall. At the wall, we outflow (copy) 
the components of velocity and magnetic field that are par- 
allel to the wall, and mirro r the perpendicular components 
dTchekhovskov et ai1l2009al) . For a given BH spin and grid 
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resolution, we choose the value of Rq such that there are 7 to 
10 grid cells between r; n and ru- This ensures that the inner 
radial grid boundary r = n n is causally disconnected from the 
region outside the BH horizon. 

At time t = we initialize the simulation with a purely 
poloidal field configuration as desc ribed in equation (or 
equation |3] in the case of the IBZ771 paraboloidal model) and 
we run the simulation until tf = max(100M,20/Oe, 10ro). Be- 
cause of the dragging of frames by the spinning BH, the mag- 
netic field develops a toroidal component which propagates 
out along field lines at nearly the speed of light. Behind this 
outgoing wave, the solution settles down to a steady state and 
we study the properties of this steady solution. We have veri- 
fied by running selected simulations for 10 times longer than 
our fiducial tf that the near-BH regions of our numerical solu- 
tions have truly reached steady state. 

3. RESULTS 

3.1. Power Output of Black Holes with Razor-thin Disks 

As explained below equation all our jet models corre- 
spond to the case of a razor-thin disk. We describe here the 
re sults w e obtain for these models. 

IBZ77I showed that the luminosity of a force-free jet from a 
slowly spinning BH (a <C 1), embedded in a regular magnetic 
field with a fixed total flux (sourced by toroidal currents in a 
razor-thin disk), is proportional to the square of the BH spin 
and the square of the magnetic field strength at the horizon. 
If we include the length scale of the horizon 2M to obtain the 
correct dimensions, we may write (the choice of the numerical 
prefactor will become clear below) 

where $ tot oc BM 2 is the total poloidal magnetic flux in the 
jet, and A: is a constant which depends on the field geometry, 
e.g., k = k mono = (6ir)~ l w 0.054 for a mon opolar field (v = 
0) and k = k pam ss 0.044 for the paraboloidal |BZ77 l ge ometry 
(equation |3J. We refer to equation (0]i as the original IBZ77I 
scaling. 

Because equation (@]i was derived in the limit a <C 1, it can 
be extended to larger values of a in several ways. In particular, 
we could replace the length scale 2M by the horizon scale r^, 
where r^ is defined in equation dTJ. In fact, this is a more 
natural scaling since the angular frequency of the BH, 

n H (a)=-^— (5) 
2r H (a) 

clearly plays an impor tant ro l e in determining the power i n 
the je t at the horizon ( lBZ77t iMcKinnev & Gammiel 120041) . 
IBZ77I found that, for a fixed fie ld strength at the BH horizon, 
the po wer in the outflow obeys (lBZ77blMcKinney & Gammiel 
12001 : 

/ , ocO(f2 H -0)| r=m ocO^ (6) 

where, based on dimensional argument, the field line angular 
frequency, J7, is propor tional to £7h (see A ppendix lAV Nu- 
merical simulations by iKomissarovl (1200 lb showed that the 
IBZ77I effect achieves maximum efficiency when fi sa flu/ 2. 
The result was found to be true for a = {0.1,0.5,0.9}, demon- 
strating that this result is valid even in the non-linear regime. 
Now, replacing the length scale 2M with the horizon scale re 
in the expression for power we obtain: 

p BZ2 (ao=£$L^H- (7) 



In Appendix [A] we derive this formula analytically from first 
principles. The scaling of jet power (0 was confirmed in the 
numerical simulations by Krasnopolsky (private communica- 
tion). In general, k is a constant factor whose value depends 
on the field geometry near the BH. For a slowly spinning BH, 
a ~ 4fluM, and equation (0 reduces to the standard IBZ77I 
scaling which was derived in the limit a <C 1. However, as 
expected from the above discussion and as we will confirm 
shortly, equation (0 is a better approximation for higher spins 
and is quite accurate up to a w 0.95, beyond which it requires 
a modest corr ection . We refer to equation (0 as the modified 
second-order IBZ771 scaling, or simply the BZ2 scaling. We 
classify the order of a scaling by the maximum power of flu 
up to which the scaling maintains its accuracy. As we will see 
below, for large spins a ~ 1, scalings higher than the 2nd or- 
der are required to obtain good agreement with the numerical 

re sults. 

iTanabe & NagatakH (120081) found 4th order corrections to 
the BH power output by performing the expansion in powers 
of BH spin a. As we have argued, a more accurate expansion 
is in powers of the BH rotati onal frequency flu (equatio n|5b. 
Recast in powers of fin, the ITanabe & Nagataki (120081) ex- 
pansion becomes: 

P BZ4 (ft H )«£$L(^ + a^), (8) 

where the 4th-order coefficient a = 8 (67-67T 2 ) /45 w 1.38. 
We analytically derive this formula from first principles in 
Appendix|B] Note that the expansion only contains even pow- 
ers of f^e since the power is independent of the sense of BH 
rotation. As we will see, this 4th order correction agrees well 
with the numerical results and is an improvement over the 2nd 
order formula (0. However, at high spins, a > 0.99, even the 
formula (J8j becomes inaccurate. Below we present a more 
accurate 6th order formula (see equation|9]l. 

We have performed numerical simulations of force-free jets 
confined by a rigid wall (as described in © for a wide range 
of field geometries and BH spins. We numerically explored all 
possible combinations of v = {0,0.25,0.5,0.75, 1}, ro = M x 
{0, 1,5, 10, 100} and a = {0.1,0.2,0.3,0.5,0.7,0.9,0.94,0.97, 
0.99,0.99769,0.999,0.999769,0.9999}. We also performed 
simula tions in which we started the field with the paraboloidal 
IBZ77I geometry d3j. 8 In addition to force-free simulations, 
which neglect plasma inertia, we have also performed MHD 
(mass-loaded) simulations for selected field geometries and 
spins: v= {0,1}, r = M X {0, 1, 10} and a = {0.1,0.9,0.99, 

0. 9999}. 

Figures [TJ and [2] illustrate the effect of the shape parame- 
ter v and the transition radius ro on the jet geometry. The 
larger the value of v, the more collimated is the jet. The larger 
the value of ro, the more monopolar-like is the jet geometry 
near the BH. The shape of the poloidal field lines weakly de- 
pends on the BH spin: careful examination of the figures re- 
veals that the field lines tend to come closer to the jet axis 
for faster spins, and this effect is largest near the BH hori- 
zon. The magnetosphere clearly divides into an outflow re- 
gion (u r > 0) and an inflow region (u r < 0) separated by a 
stagnation surface at which the radial velocity vanishes. This 

8 We note that this field geometry and the paraboloidal geometry with v = 

1 , ro = are inherently difficult to study numerically: the wall in these models 
makes such a small angle with the surface of the BH horizon that there exists 
no physical solution for the velocity in the immediate vicinity of the w all. We 
have obtained numerical solutions corresponding to the paraboloidal IBZ77I 
model only for a < 0.9 and the paraboloidal v = 1 model only for ro > 1. 
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is similar to what is seen in sim ulations of magnetized turbu- 
lent tori around spinning BHs dMcKinnev & Gammid 120041 : 
!McKinnevl2006bt iMcKinnev & Blandfordll2009l) . 

Figure |3 shows the numerically measured power output of 
all our models as a function of the BH horizon frequency H 
(defined in eq[5]>. The colored stripes correspond to the three 
scalings: BZ (equation[4]i, BZ2 (equation[7]l, and BZ6 (equa- 
tion [9]). For a given BH spin, a monopolar field geometry 
(v = 0) produces a more powerful jet since it has a larger value 
of the pre-factor k mono rj 0.054 compared to the paraboloidal 
IBZ77I geometry (equation[3 fc para w 0.044. IBZ77I) . This differ- 
ence in k determines the width of the colored stripes in Fig. [3] 
The m odel given in equation (O with v = 1 is close to the 
IBZ77I paraboloidal solution and has nearly the same power. 
Models with intermediate values of v (0 < v < 1) or with non- 
zero values of ro have power output intermediate between the 
two limiting models and form the vertical clusters of numer- 
ical points at each a in Figure [3] Independent of the value of 
v, we find that jets with ro much larger than the outer "light 
cylinder" 9 radius ~ 1 / Oh have luminosities very similar to 
that of a monopolar jet. 10 

As e xpected, for small BH spins, a < 0.3, both the original 
IBZ77I scaling (equation |4]i and the second order BZ2 scaling 
(equation 17) agree well with the numerical results. As we 
go to higher spins, the BZ2 scaling (O continues to follow 
the simulation data points accurately, while the original IBZ77I 
scaling under-predicts the power (e.g., by a factor « 3 at 
a = 0.99). 

Careful examination of Fig. [3^ and especially of the inset 
Fig.|3ja, which shows a blowup of the a — > 1 region of the plot, 
reveals a flattening in the numerical data points at spin values 
a > 0.95: the BZ2 formula (0 over-predicts the jet luminosity 
by about 25% as a — > 1 (in agreement with R. Krasnopolsky, 
private communication). It is not clear that this corner of pa- 
rameter space is particularly relevant for astrophysics, nor is 
the effect very large. Nevertheless, for completeness we note 
that the flattening of the jet power can be well-modeled by in- 
cluding higher-order corrections to the BZ2 formula (O, e.g., 
by the BZ6 formula which we derive in Appendix|C] We give 
here a simplified version of this formula: 

^ Z6 (0 H )«*C(0 H + aO H +/30 H ), (9) 

where the value of a is determined analytically, a « 1.38 
(same as in the BZ4 expansion, equation|8]l and /3 is found nu- 
merically by least-square fitting equation © to the full BZ6 
analytic formula derived in Appendix [C] ,8 « -9.2. The gray 
stripe in Figure [3] compares this formula to our numerical re- 
sults for the full range of models. Higher order corrections 
have no effect at small spin values (the gray and light red 
stripes lie on top of each other) but do a good job of reproduc- 
ing the flattening in the jet luminosity at spin values a > 0.95 
and the slight but systematic increase in the power output of 
the numerical jets above the light red stripe at a < 0.9 (this 
increase is especially apparent in Figure [3J)). In anticipation 
of future discussion, it is useful to express the power at low 

' By the "light cylinder" we mean the Alfven surface. 

10 This highlights the fact that jet power output is set by the field line shape 
close to the BH, i.e., inside the light cylinder. It suggests that communication 
along the jet is maintained by Alfven waves (rather than fast waves), so that 
the outer light cylinder, which acts as a sonic surface for Alfven waves, pre- 
vents signals propagating back to the BH from further out. As a result, the 
shape of the wall or the properties of the confining medium outside the light 
cylinder have no influence on the power output in the jet. 



spin in terms of the maximum achievable power at a = 1 : 

P(a)~0.32fl 2 P(fl= 1), a<0.3. (10) 

We now look into the origin of the differences in the power 
outputs of the various model jets, as well as of the numeri- 
cal trends discussed above. We focus on two limiting cases: 
monopolar jet (v = 0) and paraboloidal jet (v = 1, ro = 1). 

First, let us recast the power output of the jet in a conve- 
nient form. In a stationary axisymmetric force-free flow, sev- 
eral quantities are conserved along poloidal field lines (de- 
fined by $ = const). Two of these are the field line an- 
gular velocity O('I') and th e enclosed poloidal current I(^>) 
dTchekhovskov et al. 2008). The power output of a force- 
free jet may be written as the integral of the outward Poynt- 
ing flux S r = -QB r B v over a spherical jet cross-section 11 . In 
this notation, the lower component of the toroidal magnetic 
field is up to a numerical factor the enclosed poloidal current, 
-2irB v = Using this notation, which is very similar in 

appearance and meaning to the usual special relativistic nota- 
tion, we obtain the total power outp ut of the BH by integrating 
over the surface of the BH (see also lBZ77h : 

P = J J S r dA = 2^ ' WldA^lj Q(^!)I(^)d^, (11) 

where the field strength B r times the area element dA gives 
the magnetic flux through that area, d$> = 2-nd^ = B r dA, and 
the numerical factor of 2 accounts for the two hemispheres of 
the BH. 

Figure @] shows for a monopolar jet the angular profiles 
of angular velocity uj = O/Oh, enclosed poloidal current i = 
//Oh, and power output p = P/0 H . The particular scalings by 
Oh have been selected based on equation (Q so as to remove 
any obvious trends as a function of spin. This allows us to 
compare models with different spins on the same scale. At 
low spin, a < 0.1, we have excellent agreement betw een the 
numerical models and the analytic solution obtained bv lBZ77l 
shown by the dotted lines. For larger spins up to a < 0.9, 
both the dimensionless angular field line velocity lj and the 
enclosed current i increase with increasing a (Figures @^,b). 
According to equation (fTTI) . this should result in an increase 
in the normalized jet power p, as confirmed in Figured};. This 
is the reason for the small but systematic increase in jet power 
above the estimate Q. For a > 0.95, we find that lj, i, and 
p all decrease relative to (0, causing a flattening of the jet 
power at these extreme spins. The reason for the decreased 
power is related to a change in th e pol oidal field geometry of 
the jet near the BH horizon (see 33.21 and Appendices IB1 and 
let : while at low spins the magnetic field is nearly uniform 
across the jet for all of our models, at high spins the poloidal 
field becomes non-uniform with a maximum field strength at 
the jet axis and a minimum near the wall. Since it is the field 
geometry near the BH that sets the power output (see foot- 
note [TO] and Appendix IaI). it is logical that these changes in 
the field geometry l ead to changes in the power output. We 
demonstrate this in j33.2l (see also Appendices IB1 and [CT> . 

1 1 Here the GR notation is simplified and appears like the non-GR expres- 
sions (apart from some sign conventi ons) by using the notational conventions 
in appendix B of iMcKinnevI )2005aT) and in lMcKinneyl )2006aD . In this no- 
tation, B l = *F" , Bj = "Fu, Ej = F u /y/=g, E' = F"^, and fi = -E e /B r , 

where F is the faraday tensor, F is the Maxwell tensor, and \/—g is the square 
root of minus the determinant of the metric. Horizon surface area elements 
are given by dA = y^gdddcji. 
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Figure 4. Angular dependence of various quantities in a monopolar jet (u = 
0) as a function of the poloidal flux function The different curves in each 
panel correspond to different values of the BH spin (see legend). From top 
to bottom the panels show the normalized field angular velocity uj = f2/f2H, 
the normalized enclosed poloidal current i = I/Qn, and the normalized jet 
luminosity p = uii = QI/Q^ = P/H^. The analytic TBZ77I solution, shown 
with dotted lines, provides an excellent description of the numerical results 
for all spin values a < 0.99. Beyond this value of a, the quantities ui, i, 
and p all become noticeably smaller than the analytic solution. This trend is 
removed by the BZ6 solution (9), as shown in FigurefJ] 

Figure [5] shows that paraboloidal jets exhibit very similar 
trends with increasing spin as their monopolar counterparts. 
The differences are in details, e.g., the angular velocity pro- 
file (Figure[5j:) is now no n-uniform even for low spin values, 
as predicted by the lBZ77l analytic solution shown in the figure 
with dotted lines. The agreement with the analytic solution is 
not as perfect as for the monopolar model, but this is because 
the poloidal field line shape near the BH in our v = 1 , ro = 1 
paraboloidal jet differs sli ghtly f rom the IBZ77I paraboloidal 
shape. For our numerical IBZ77I paraboloidal jets the agree- 
ment with the analytic solution is very good. 

In all our numerical jets the conserved quantities and 
r2($) are preserved along field lines to better than 10%. We 
reran a selection of models at twice the fiducial resolution in 
both the radial and angular directions. We found differences 
of less than 5% in the total power output, indicating that our 
models are well-converged. 

3.2. Power Output of Black Holes with Thick Disks 

In all the models we described so far, the base of each polar 
jet covered a full 2ir steradians at the BH horizon. However, 
observational evidence strongly suggests that low-luminosity 



Figure 5. Similar to Figure|4]but for a paraboloidal jet (v = 1, tq = 1). Com- 
parison with Figure[4]shows that, for the same BH spin, the angular velocity 
uj is smaller and the enclosed current i larger in a paraboloidal jet compared 
to a monopolar jet. These two effects combine to give a smaller power out- 
put p = toi in the paraboloidal solution. Note that, whereas a monopolar 
jet rotates more or less like a rigid body, a paraboloidal jet has a significant 
variation of u) across its cross-section. 

BHs (A < .01, {Q} are surrounded by thi ck accretion disks 
or ADAFs (Naraya n & McClintockl 120081) with thicknesses 
H/R ~ 1. Here and below by the "disk thickness" H/R we 
mean the angular extent at the BH horizon of the region exte- 
rior to the Poynting-dominated jet, i.e., the total thickness of 
the gaseous disk plus any magnetized corona or heavily mass- 
loaded wind above the disk. When a BH is surrounded by a 
thick disk/corona, equatorial field lines from the BH at lower 
latitudes pass through the disk/corona, become turbulent and 
produce a slow baryon-rich wind, whereas polar field lines 
at higher latitudes lie away from th e disk gas and pro duce a 
Poynting-dominated relativistic jet dMcKinnevll2005bl) . How 
does this effect modify the dependence of jet power on BH 
spin? 

Assuming that both the total magnetic flux $ to t threading 
the BH horizon and the angular thickness H/R of the ac- 
cretion disk/corona are independent of the BH spin, we can 
compute the power that is emitted in the Poynting-dominated 
region of the jet. We assume that the models we have de- 
scribed earlier continue to be valid, except that we integrate 
the jet power only over field lines that cross the horizon out- 
side the ±H/R zone of the disk/corona. This procedure is 
well-motivated by general relativistic magnetohydrodynami- 
cal simulations of thick accretion disks which show that the 
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Figure 6. Jet power output of various models as a function of dimensionless BH horizon spin frequency i^h = ^h/^hO? =1) (see eq.|5J for four different choices 
of disk thickness: razor-thin disk with H/R = 0, thicker disks with H/R = 0.5, 1, and 1.25. For reference, the following values of BH spin, a = {0.2,0.5,0.9, Ik 
correspond to approximately the following values of wh: {0.1,0.27,0.63, 1}. Numerical results for different models are shown with different symbols (see Fig.[3] 
caption for details). The second order analytic BZ2 solution is shown with dotted lines, the 4th order BZ4 solution with dashed lines, and the 6th order BZ6 
solution with solid lines. As the disk becomes thicker, the spin dependence of jet power becomes steeper. This steepening is most accurately described by the 
BZ6 formula. Note that this formula also reproduces the flattening of the jet power for razor-thin disks (H/R = 0) at high spins (Figure|5J. For reference, we also 
show the slopes of various power-law scalings P <x Q'^ with straight line segments. 



relativistic jet subtends a well-defined solid angle for a given 
gas pressure scale height dMcKinnev & Gammid 120041) . We 
consider models with H/R = 0.5, 1, and 1.25, and compare 
the results with the case of a razor- thin disk (H/R = Q). 

The results for the jet power as a function of disk thickness 
and spin are shown in Figure [6] As we have already seen, 
the jet power for a razor-thin disk scales as P oc (this scal- 
ing is shown with dotted lines) until a < 0.95 after which it 
levels off. As the disk becomes thicker, the scaling changes 
qualitatively. For all thicknesses, H/R = 0.5, 1, and 1.25, the 
power dependence on the spin follows the same Qj, power- 
law at low spins. However, at higher values of a, the power 
increases more steeply. The break occurs roughly at a ~ 0.7, 
with a moderate dependence on the disk thickness. Above 
the break for the case H/R= 1 we have P oc fl^ and for the 



case H/R= 1.25 we have P oc f2„. This steep depe ndence is 
similar to what was observed by [McKinnevI (12005b!) in his nu- 
merical simulations of turbulent accreting tori. 12 We observe 
the same steeppower dependence also in our ideal GRMHD 
simulations (H). 

We now explain the reason for the steep dependence of 
jet power on flu as a — > 1 . We saw in Figures Q] and [2] 
that as the BH spin increases, magnetic field lines rearrange 
laterally and concentrate aroun d the axis of rotation (see 
iKomissarov & McKinneyll2007l for an explanation of this ef- 
fect in terms of hoop stresses). Figure [7] shows that this 
leads to a non-uniform distribution of radial magnetic field 
B r across the jet, with B r having a maximum near the rota- 

12 Note that only even powers can enter the expansion of the jet power in 
terms of Qu because the power is an even function of £7h- 
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tion axis and a minimum near the jet boundary. Since the 
electromagnetic energy flux of a BH i s proportional to (B r ) 2 
at the BH horizon (see equation lAU . the concentration of 
magnetic flux near the rotation axis leads to a progressively 
larger fraction of the total energy output of the BH to be 
emitted in the polar region, giving a steeper dependence of 
jet power on the BH spin in the presence of a thick disk 
(H/R ~ 1). In a related conte xt, iMacDonakH (11984 and 
iKomissarov & McKinnev (2007|) have studied how magnetic 
flux is pulled in by a spinning B H by considering magnetic 
hoop stresses. [MacDonaldl (119841 ) appears to have missed the 
strength of this effect by mostly investigating models with rel- 
atively small a < 0.7 and by primarily looking for a change in 
the total magnetic flux accumulated at the horizon rather than 
measuring the flux redistribution on the horizon. 

We now describe an analytic approach for understanding 
the steep dependence of jet power on BH spin for thick disks 
with H/R ~ 1. While the split-monopolar magnetic field is 
an exact solution for non-spinning BHs, for spinning BHs 
the dragging of frames induces a spin-dependent perturba- 
tion to the split-monopolar magnetic field geometry. It is 
this perturbation that causes field lines to move preferentially 
toward the r otation al axis. By performing an expansion in 
powers of a, IBZ77I determined this perturbation of the mag- 
netic fiel d geometry to the lowest (sec ond ) order in BH spin a 
(see a lso lMcKinnev & Gammiel2004l and lTanabe & Nagatakil 
|2008|). Accounting for this p erturbation in the field geome- 
try, tfanabe & Nagatakil (120081) deter mined th e BH power out- 
put more accurately than the original IBZ77I derivation, to the 
4th order in BH spin a. As we noted in § 13.11 expansions in 
powers of BH frequency flu are more accurate at high spins 
than expansions in powers of the BH spin a. Therefore, in 
Appendix [B] we perform an equivalent expansion in terms of 
Oh- 13 Figure [6] shows with dashed lines the power output 
of our ^H-based model, which we refer to as BZ4. Clearly, 
it provides a more accurate approximation for the power of 
our numerical jets than the 2nd order accurate BZ2 solu tion 
shown with the dotted lines. However, as we saw in §13. H and 
as is clear also from Figure [6] the BZ4 solution still does not 
capture a few important effects: (i) for BHs with razor- thin 
disks, it does not capture the flattening of the power output at 
a > 0.95 and thereby over-predicts the numerical results by 
as much as 25%, and (ii) for thick disks this solution under- 
predicts the numerical power by as much as 70%. 

To improve our analytic model, we have used the results 
of our numerical simulations to determine higher-order cor- 
rections to both the field geometry and the power of the BH 
at high spins. Firstly, we have obtained a higher-order accu- 
rate numerically-motivated expansion in powers of 51h for the 
magnetic field B r at the BH horizon (Appendix [C] provides 
the details). Shown with dotted lines in Figure [JJ for a wide 
range of a, this analytic approximation for B r accurately re- 
produces the numerical angular profile of magnetic field for a 
wide range of polar angles, 6*h > 0.3 or H/R < 1 .3. Secondly, 
based on this higher-order magnetic field profile, we have an- 
alytically computed the 6th order accurate approximation for 
the BH power output, which we refer to as model BZ6 (see 
Appendix ICl for more detail). Figure [6] shows the results with 
the solid lines. Not only does the BZ6 expansion correctly 
capture the flattening of the BH power at high BH spins for 
razor-thin disks, it also provides a significantly more accurate 
approximation to the jet power output for thicker disks. For 

13 This expansion reduces to the expansion {8) in the limit of H/R—} 0. 
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Figure 7. Radial contravariant component of the magnetic field strength B' 
evaluated at the BH horizon in a monopolar model (y — 0, tq — 0) as a function 
of polar angle 9u, for different values of BH spin (see legend). As the spin 
of a BH is increased, magnetic field lines progressively bunch up toward the 
BH rotation axis, resulting in an increased magnetic field strength close to 
the axis at small 6>h (this effect can also be seen in Figs.[TJ{2)- Dotted colored 
lines show the high-order analytic solution, while the various other lines show 
the numerical solution (see legend). 



instance, for a thick disk with H/R=l .25, the BZ6 expansion 
is about a factor of 20 more accurate than the BZ4 expansion. 

4. DISCUSSION 

Figure [3] shows that, regardless of the geometry of the con- 
fining wall, the total power output of a magnetized relativis- 
tic spinning BH with a razor-thin disk varies as P = A:$ t 2 ot il^ 
(equation [7]), where $ tot is the total magnetic flux threading 
the BH horizon, £1h = a /(2rn) is the BH horizon frequency, 
and ru=M [l + (1 - a 2 ) 1 ' 2 ] is the radius of the BH horizon in 
units of G/c 2 . This modified BZ2 scaling is slightly steeper 
than the original lBZ77l scaling P = k^ 2 ot (a/4M) 2 (equation^. 
A more accurate BZ6 scaling (equation [9]) accurately repro- 
duces the power-spin dependence for all values of BH spin a, 
including the limit a — > 1 . 

In the context of the radio loud/quiet dichotomy of AGN, 
following lSikora et al.1 d2007l) let us assume that supermassive 
BHs in elliptical galaxies, which manifest themselves as radio 
loud AGN, have higher spin parameters with a median a ~ 
0.9, while the BHs in spira ls, the radio quiet AG N, have a 
lower median a ~ 0.3 (e.g.. IVolonteri et aUl2007l) . Figure [3] 
then suggests that the jet powers in the two classes of objects 
(assuming similar values of <£> t ot) would differ by a factor ~ 
20. However, radio loud AGN and radio quiet AGN differ i n 
their radio luminosities by a factor ~ 10 3 dSikora et al.ll2007l) . 
What could be the reason for such a large dichotomy in jet 
power? 

One motivation for the present study was to investigate 
whether there is any strong non-linearity in BH physics that 
might cause the jet power to increase very rapidly as the BH 
spin approaches unity. If this were the case, one could pursue 
a scenario in which radio loud AGN are associated with nearly 
extremal Kerr BHs. Unfortunately, our numerical results in- 
dicate that non-linearity hardly helps. Because the total BH 
power output scales as (equation |5]l rather than simply as 
a 2 , there is a slightly steeper increase of power with a as the 
spin approaches unity. However, the scaling actually becomes 
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shallower once a increases above ~ 0.99. We have carefully 
checked the convergence of our models and we are confident 
that the results are not affected by numerical errors. There- 
fore, there is not much room for increasing the power of radio 
loud AGN jets by pushing a arbitrarily close to unity. 

Therefore, since there is not much wiggle room at the ra- 
dio loud end, we need to postulate that radio quiet AGN 
have very low values of a, say a ~ 0.03. This is uncom- 
fortably low. It is certainly feasible for an occasional BH 
to be spinning so slowly, but to have an entire population 
of BHs (radio quiet sources in spirals) with a median a of 
order 0.03 seems far-fetched. It would require spiral galax- 
ies not to have experienced any significant mergers. Further- 
more, the BHs in their nuclei should have accreted mass en- 
tirely through minor mergers with smaller companions, each 
with a tiny mass and with a random orientation of angular 
momentum dHughes & Blandfordll2003t iGammie et alj|2004t 
iBerti & Volonterill2008l) . ~ 

The second motivation for the present study was to investi- 
gate if the geometry of the confining funnel along which the 
jet propagates, which may be different for rapidly-spinning 
and slowly-spinning BHs, could lead to a substantial change 
in the jet power. This too turns out not to be the case, within 
the context of razor- thin disks. We have tried a wide range of 
geometries for the jet, as described in ^2](see also Figs. [T||2), 
but the jet power varies by no more than 20% for a fixed BH 
spin and magnetic flux. 

Our third motivation was to in vestigate the correctness of 
the results of lMcKinneyl d2005bl) who obtained a steeper de- 
pendence in the jet power compar ed to the total BH power. 
For this purpose, we considered in jj3.2l an additional geomet- 
rical effect, viz., varying the solid angle subtended by the base 
of the jet. Such a variation is expected if the accretion disk is 
geometrically thick. As an example, let us consider the re- 
sults corresponding to a BH surrounded by a disk of angular 
extent H/R=l, and let us further assume that radio loud AGN 
have spin parameters very close to unity and radio quiet AGN 
have more modest values of a. Then Figure [6] shows that it 
is possible to explain the radio loud/quiet dichotomy, i.e., a 
factor of 10 3 in jet power, if radio loud AGN have a — > 1 and 
radio quiet systems have a ~ 0.15. This is a lot more com- 
fortable than the requirement a ~ 0.03 that we found earlier 
for a razor-thin disk. Indeed, if the disk/corona is even thicker 
than H/R ~ 1, which is not unreasonable 14 , the effect is even 
stronger (see Figure |6}, and the radio loud/quiet dichotomy 
can be explained with quite modest changes in spin. 

There are several other effects that we did not con- 
sider which might either enhance or diminish the effect of 
rapid spin. For instance, we assumed that the total mag- 
netic flux threading the horizon is a constant, indepen- 
dent of BH spin. Any mechanism that enhances the to- 
tal magnetic flux threa ding the horizon would lead to a 
larger BH power output dMcKinnev 2005b; Hawlev & Krolikl 
20061 IKomissarov & McKinnevI 120071; Reynolds et al.1 120061 



Garofaloll2009l) . For example. iMcKinnevI d2005bl) found that 
the magnetic flux across the entire horizon scales as <& tot cx 

1 /2 

J7 H , so the total power will increase by another factor of 
&t ot oc fin, which further diminishes the range of BH spins 
required to explain the radio loud/quiet dichotomy. At first 

14 We note that H/R in this context refers to all the gas-dominated regions 
of the flow: the accretion disk proper, the corona and the disk wind. The net 
half-angle subtended by all these components could equal a radian or more 
in the case of a thick accretion flow. 



sight, it would appear that $ tot is determined by conditions 
far from the BH, e.g., the net magnetic flux of the gas sup- 
plied to the accretion disk on the outside, and the ability of 
the disk gas to transport this field in. However, once the field 
has been transported to the center, two general relativistic ef- 
fects kick in. First, the spin of the BH determines the size of 
the plunging re gion, and larger plunging regions can accumu- 
late more flux dGarof alo 2009). Note that magnetic flux can 
also be transpor ted to the BH through the corona outside the 
accre tion disk dRothstein & Lovelace! 120081: iBeckwith et al.l 
120091) . Second, the frame-dragging of space-time near the 
BH within the ergospheric region leads to currents that gen- 
erate hoop stresses pulling magn etic flux toward the horizon 
dKomissarov & McKinney||2007l) . These two effects are non- 
trivially coupled, although for prograde BH s pins the hoop 
stresses appear to dominate (lMcKinnevl2005bl) . while for ret- 
rogra de BH spins the size of the plunging region may dom- 
inate dGarof aloll2009l) . Another possibility is that a stronger 
dependence on spin could occur if some field lines attach be- 
tween the disk and the BH dYe & Wang||2005l) . although such 
configurations are not seen in general relativis tic magneto- 
hydrodynamical simulat ions of accretion disks dHirose et al.l 
l2004tlMcKlrinevll2005bl) . 

Yet another possibility is that mass-loading of AGN jets 
might have a large effect on the jet power. For exam- 
ple, the IBZ77I mechanism only operates at sufficiently high 
magnetization for a given BH spin. Thus, the magne- 
tization and BH spin can together introduce a "magnetic 
switch" mec hanism that can trigger powerful jet formation 
from t he BH dTakahashi et alfl990t iMeier et alj[l997t iMeierl 



119991; IKomissarov & Barkovl 20091) . Another type of mag- 
netic switch can be due to changes in the field geometry 
from dipolar to multipolar, wh ich leads to significant mass- 
loading of the polar regions (jMcKinnev & Gammid 120041 ; 
IBeckwith etail 120081 iMcKinnev & Blandfordl 120091) and a 
factor of ~ 10 weaker total BH power output. The physics 
of jet mass-loading is presently uncertain, and this question 
needs to be investigated in more detail. Finally, changes in 
the disk thickness may also result in differences in the amount 
of magn etic flux accu mulated or generated by turbulence near 
the BH dMeierll2001l) . 

Finally, we note that we have implicitly assumed in this 
paper that the radio luminosity of a jet is directly proportional 
to the total energy flux (Poynting and kinetic power) carried 
by the jet. Perhaps this is not the case. Any non-linearity in 
the mapping between radiative luminosity and jet power could 
have important consequences. In particular, we note that the 
interstellar medium (ISM) in a typical elliptical galaxy is very 
different from that in a typical spiral. Since the radio emission 
in a jet is produced when the jet interacts with the external 
ISM, this difference may well lead to a strong effect on the 
radio loudness of the jet. 

In application to gamma-ray bursts (GRBs) and collapsars, 
an interesting question is whether they are powered by the 
BZ mechanism through an outflow fr om a central BH or by 
an out flow from an accretion disk. IKomissarov & Barkovl 
d2009l) suggest that the BZ effect i s operating in such a sce- 
nario (however, see|Nagataki 2 0091) . 

5. CONCLUSIONS 

We set out in this paper to explain the radio loud/quiet di- 
chotomy of AGN in the context of the BH spin paradigm. For 
razor-thin disks, we found that BH spin alone is insufficient to 
explain the observations even if the radio loud and radio quiet 
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populations have very different merger and accretion histo- 
ries. However, we found that the presence of a thick disk, such 
as an ADAF, can significantly enhance the spin dependence of 
the power output, to the extent that it can reasonably account 
for the observed radio loud/quiet dichoto my. Our only mod- 
ification to the revised spin paradigm of ISikora et al.l (120071) 
is that both populations should contain a BH surrounded by a 
thick disk such that the jet subtends a small solid angle around 
the polar axis. 

These results were obtained by performing general rela- 
tivistic numerical simulations of collimated force-free and 
MHD jets from spinning magnetized BHs for a wide range 
of spins (up to a = 0.9999) and jet confinement geometries. 
We showed that, regardless of the geometry, a BH threaded 
with a magnetic flux $ tot and surrounded by a razor-thin disk 
produces a jet with power P « £$ t 2 0t r^, where k is a known 
constant factor which depends only weakly on the field ge- 
ometry, re =M[1 +(1 -a 2 ) 1 / 2 ] is the radius of the BH horizon, 
and ilu = a /(2rn) is the angular frequency of the BH. This 
result gives a somewhat steeper dependence of jet power on 
a compared to the original scaling P oc a 2 obtained by |BZ77l 
Nevertheless, we conclude that for a fixed magnetic flux $ tot , 
even this revised scaling is much too shallow to explain the 
radio loud/quiet dichotomy of AGN. Our goal, therefore, was 
to identify any other effect that may cause the jet power to 
depend more steeply on BH spin. 

We found that such an effect naturally exists. We showed 



that the power output of a BH surrounded by a thick accre- 
tion disk with H/R ~ 1 (this is the effective thickness of 
the disk, corona and mass-loaded disk wind), as expected in 
systems with advection-domi nated accretion flows (ADAFs, 
iNarayan & McClintockll2008l) . is P oc Q^, and even oc Cl^ for 
very thick disks ( jj3.21 i. In this case we can explain the ra- 
dio loud/quiet dichotomy by having two different populations 
of galaxies with modestly different BH spins. For the case 
H/R= 1, the radio loud population needs to have large spins 
a ~ 1 while the radio quiet AGN population needs to have 
a ~ 0.15. Such spin values may plausibly result from differ- 
ences in the merger and accretion histories of supermassive 
BHs in elliptical and spiral galaxies (® . 

We worked out in the Appendices a first principles analytic 
model which accurately reproduces our numerical results for 
the jet power over a wide range of BH spin and disk thickness 
(Figures|3]|6]l. 
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APPENDIX 



A. SECOND-ORDER-ACCURATE EXPANSION OF BLACK HOLE POWER 

In this section we present a compact derivation of the BZ77I effect. We determine the power output of a spinning BH embedded 
into an externally-imposed split-monopolar magnetic field. This magnetic field is given by the flux function (O with v = and 
ro = 0. The main difference of this derivation is that we perform it in the powers of the "natural" variable - the BH angular 
frequency flu that plays an important role in determining the BH power output. The power output density evaluated at the 
horizon of a spinning BH is dBlandford & Znajekl[T977HMcKinney & Gamm ie 2004) 



F E (9) = [2(B r ) 2 n(fl H - 0)rM sin 2 9] 



(Al) 



where the quantity O is the angular frequency of magnetic field lines at the BH horizon and B r is the radial field strength at the 
horizon. 

In order to determine the power output (lAlt . we need to know two quantities as functions of polar angle at the BH horizon: the 
radial magnetic field, B r , and the field line angular frequency, f2. The element of the magnetic flux d& through the BH horizon is 
related to the radial magnetic field at the horizon through the following differential, 

d$ = lirdm = 2nB r ^/=gd9, (A2) 

where g is the determinant of the Kerr metric in the Boy er-Lindquist coordinates, y/—g = (r 2 + a 2 cos 2 9) sin 9. This formula closely 
resembles its cousin in the spherical polar coordinates and flat space where one has ^/-g = r 2 sin 9. More generally, in place of 
B r we can use any vector field (e.g., energy flux), and the result is the differential of that flux. Assuming that the perturbations 
to the magnetic field away from a perfect split- mono pole are higher order, we neglect them and obtain, by differencing the flux 
function <J2J (with v = and ro = 0) according to dA2b . an identical result to that in flat space: 

B r =% ot /r 2 , (A3) 

where we neglected terms of order and higher. While the distribution of O needs to be self-consistently determined by solving 
the non-linear equations describing the balance of electromagnetic fields (we do so numerically in ^3j, here we make a simple 
yet accurate estimate based on the ener gy argument. Let us assume that the system chooses such a distribution of £1 that it causes 
an extremum in BH power output ( lAlb - Such a value is clearly 

fl = fl n /2 (A4) 

since it maximizes the BH power output ( IAU (see, e.g. jBeskin & Kuznetsoval2000l) . Despite the simplicity of this estimate, it is 
remarkably close to the true solution for the split-monopolar geometry as obtained from the numerical simulations (® . Plugging 
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equations (I A3t and (IA41 > into the power output density ( lAlb . we obtain 

F E (9) = 2 (% ot r- 2 ) 2 (n H /2) 2 rMsm 2 e. (A5) 

Integrating up this power output density in angle in the same way as we integrated the magnetic field in equation dA2l ). we obtain 
the full power output into jets with an opening angle 9y. 



P = 2x2it F E (9)^g 



d$, (A6) 



where the extra factor of 2 accounts for the fact that there are two jets, one in the northern and one in the southern hemisphere. 
Note that we are inte rested in an expansion of power up to 2nd order in D,u- S ince the factor Fe(0) oc is already second order 
in fie (equation IA5b . we can without loss of accuracy evaluate the formula (IA6t at r = ru(a = 0) = 2M and replace ^/-g with 
[(2M) 2 sinfl] . After plugging into <|A6) for F E (9) using JA5b and evaluating the integral out to 9j = tt/2, i.e., computing the full 
power output of the BH, we get: 

,ir/2 

P=ir^i t n 2 1 sin 3 0d0 = 27T* 2 ot Q 2 J /3, (A7) 
Jo 

which is accurate to second order in 51 H . In terms of the magnetic flux $ tot = Itt^^, the formula becomes 

P = k®lPl, (A8) 

with k = 1/ (67r), which reproduces Q. 

B. FOURTH-ORDER-ACCURATE EXPANSION OF BLACK HOLE POWER 

In Appendix|A]we have derived a second order accurate expression for power output of the BH in terms of the hole frequency 
fin- The BH was embedded with a split-monopolar magnetic field. Let us now improve the accuracy of the previous derivation, 
this time retaining the higher order terms, up to fi^. A similar derivation was performed by iTanabe & Nagatakil (12008) but in 
powers of BH spin a. Here we derive the expansion in terms of the natural variable Sle which allows to use the expansion for 
nearly maximally spinning BHs. We also explicitly present the angular dependence of the BH power output. 

In order to obtain a hi gher-or der a pproximation to the power, this time we need to keep some of the higher terms we neglected 
in equations for B r (IA31 > and Q, (IA4b . Since Q, is an odd function of BH frequency, it contains only odd powers of fin, therefore a 
higher order approximation for it has the following form: 

n = 1 / 2 H +O(0| [ ) ) (Bl) 

where 0(17^) denotes any third order or higher order terms in £7h- Since P oc f2j^, these higher-order terms contribute to the 
power output only terms of order higher tha n Q(flt ,), therefore we neglect them. We do need, however, to include the terms that 
come from the higher order expansion of B r . IBZ77I showed that the dragging of frames by the spinning BH perturbs the magnetic 
field away from an exact split-monopole and have derived a second order correction to the flux function, \I/ r 2 , in powers of BH 
spin so that the full flux function has the form 

#0, 9) = V (9) + a 2 V 2 (r, 0) + 0(a 4 ) (B2) 
= * O (0) + 16^2M) + O(f£), (B3) 
where we have used a = 4S7h + 0(^h)- Here the zeroth and second order perturbations to the flux function are 

*o(#)= l-cos0, *20, 9) = f{r) sin 2 9 cos 9, (B4) 

where f(r) is a known function of radius r, but for further discussion only its value at the horizon of a non-spinning BH, 
f(r = 2) = (56-37T 2 ) /45, is needed (this is because in the expansion ( IB3b we formally evaluate the coefficients at flu = 0, 
r=r H = 2). 

Combining expressions ( lB31 l. (IB4b and (1A2K we obtain the 2nd-order-accurate radial magnetic field at the BH horizon: 

nr (l+4ft H 2 ) 2 [9 + f>H 2 (-49 + 6^ 2 )(l + 3cos20)] 

B = 7 : \ ■ ( B5 ) 



' (V 2 (1 +4f2 H 2 ) 2 + 16f2 H 2 cos 9 1 



Combining this expression with (lAlb . (IA4I) and plugging into dA6b . we numerically obtain the angle-dependent enclosed power 
P BZ4 (9,n H ) shown in Figure|6]as P BZ4 (9 = tt/2-H/R, f2 H ) with dashed lines. This result, expanded to 4th order in powers of 
S^h, is: 

P BZ4 (8, fn ttQh [Vi sin 4 (6» /2)(cos 9 + 2)] 

+ ttQ£ [90 (3tt 2 - 32) cos 9 + (970 - 105tt 2 ) cos 39 + 9 (3tt 2 - 26) cos 59 + 32 (67 - 6tt 2 ) ] /270+ 0(1^) (B6) 
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Clearly, at low spins, the second-order piece dominates, which we show in Figure [6] with dotted lines. However, at high spins, 
the fou rth order piece can become dominant, which is confirmed in Figure[6] To see this more clearly, we perform an expansion 
of ( IB6b in powers of disk/corona thickness, H/R = it/2-9, and obtain to second order in H/R: 

P BZ4 (H /R)kQ 2 h {2. 09 -3. UH/R + 0[(H/R) 3 ]} + n 4 H {2. 9 +1.7H/R + 0[(H/R) i ]}+0(^), (B7) 

where for clarity we have numerically evaluated the coefficients to two decimal places. This expansion makes it clear that as H/R 
increases, the relative importance of the fourth order term increases. This also explains why in Figure[6]the P oc £1 H dependence 
becomes more prominent for larger values of H/R as opposed to smaller values. 
Finally, we note that at the midplane the power output takes the following form (expanded up to 4th order in fin)- 

p BZ \e = tt/2) « 2 7 3 *r , [n 2 H +n%8 (ei-en 2 ) /45] +o(0 H ) = 2 y 3 * 2 ol [n 2 H+a n 4 H ] , (B8) 

where a = 8 (67-6-7T 2 ) /45 « 1.38. In this formula we have reintroduced \l/ tot which was set to unity for the rest of this section. 
This result can also be expressed in terms of the total flux in the jet using $ tot = 27r5 , t ot- 

C. SIXTH-ORDER-ACCURATE EXPANSION OF BLACK HOLE POWER 

Figure [6] shows that the fourth-order BZ4 solution for power is more accurate than the second-order solution. However, at high 
BH spin, a > 0.95, it requires a more than a factor of 3 correction in order to reproduce the numerical solution. Also, the fourth 
order BZ4 solution does not reproduce the flattening of the power dependence on the BH spin for razor-thin disks (H/R = Q) at 
a > 0.95. 

Inspired by the success of the previous section, we would like to derive a sixth-order-accurate expression for the power. 
However, for that we would need to know the expansion of the flux function to the fourth order and of the field angular frequency 
to the third order. None of these are known analytically, therefore, we adopt a numerical approach. 

Figure^ shows the angular profiles of the field line rotation frequency O for different BH spins. While the deviations from the 
zeroth order approximation O = Oh/2 are present, their relative magnitude is very small, < 10%. These 10% changes in O/Oh 
translate into at most 1 % changes in the power output because the power depends quadratically on the magnitude of the higher 
order correction (0- Oh/2) (see equation[6]). We are interested in the corrections of order ~ 10-70% (the level of inaccuracy of 
the BZ4 solution), therefore we neglect the higher-order corrections to Oh in deriving the sixth order solution. 

The corrections to the magnetic field shape are, however, dramatic. Figure [7] shows the angular distribution of the radial 
magnetic field B 1 ' at the BH horizon. As the BH spin increases, B' develops a progressively large non-uniformity in angle and 
deviates from the 4th order solution by factors of a few. We therefore, attempt to find a numerical fit to the angular magnetic field 
dependence at the BH horizon for a = 0.9999 by fitting to it the following trial function: 

# = *o + 160 H * 2 + H * 4 , ( C1 ) 

where the first two terms are given by equations dB4l ). We look for the spin-independent part of the third term, ^4, in the following 
form 

* 4 (6') = sin 2 (6l)[ciCos Ql 6l + C2Cos Q2 6i + C3Cos Q3 6l + C4Cos Q4 6l] , (C2) 

where we choose a\ = 25, oli = 7, 03 = 3, 0:4 = 1. In order to determine 4 coefficients C1-C4, we match the numerical solution for 
B r at a = 0.9999, shown in Figure[p at 4 angles: 9 = 0.1, 0.5, 0.7, tt/2. We find c { « 26.16, c 2 « 22.72, c 3 « 13.54, c 4 « 2.08. 
Figure [7] shows as dotted colored lines the solutions due to jC 1 b and dC21 i, with the above values of expansion coefficients, for 
various values of BH spin. Clearly, the analytic fit is a very good match to the power output at 9u > 0.3. Very close to the rotation 
axis, however, (at angles smaller than 0.3) the fourth-order-accurate solution to the flux function jC 1 b is not enough: while we 
have a nearly perfect match between our fit and the numeric solution for B6 at high (a = 0.9999) and mid-range (a = 0.5) spins, 
our fit to B' deviates by up to ~ 25% at the in-between spins (a ~ 0.9-0.99). However, the total power emitted into the range of 
polar angles 9 < 0.3 is very small, therefore these deviations of our fit from the numerical solution hardly influence the fit to the 
power output. 

We also considered a dir ect fit to the vector spherical h armonic functions that form a complete set for the vector potential 
as given by equation B8 in iMcKinnev & Naravanl d2007al) . However, even an expansion up to I = 10 did not fit the very steep 
behavior of B r near the polar axis. However, otherwise, even only using up to I = 6 does a reasonable job at fitting the numerical 
solution for the total power vs. spin and angle. This fact and the fact that a power of 25 for cos 6* was required to fit the numerical 
results demonstrates that the numerical solution at a ~ 1 is highly non-linear with respect to 9 and would be quite difficult to 
derive analytically. This proves the usefulness of the numerical simulations. 

Now we are in a position to analytically compute the high -ord er-accurate power of our jets. Using dA2l ). we compute the r adial 
field on the h orizo n, B r , from the fourth-order flux function jCl b . We then insert this field and the field angular frequency O dA4b 
into formula iA6i and obtain the angular-dependent jet power. This power, which we refer to as the sixth-order analytic BZ6 
solution, is shown in Figure [6] with solid lines for various disk thicknesses and spins. (We compute these lines by numerically 
integrating the analytically-determined power in our jets. We note while formally this solution is 6th order-accurate, its expansion 
in powers of Oh contains important terms up to 10th order. This highlights the non-linearity of the problem.) Clearly, these lines 
approximate the numerical data points very well, within 20% for the whole range of disk thicknesses and spins that we have 
explored. Over most of the parameter space the errors are smaller than this value (they are largest for the thicker disks with 
H/R>\ .25 that have the BH spin in the range 0.8 < a < 0.99). 
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